Monte Carlo studies of extensions of the Blume-Emery-Grifnths model 
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We extend the Blume-Emery-Griffiths (BEG) model to a two-component BEG model in order to 
study 2D systems with two order parameters, such as magnetic superconductors or two-component 
Bose-Einstein condensates. The model is investigated using Monte Carlo simulations, and the 
temperature-concentration phase diagram is determined in the presence and absence of an external 
magnetic field. This model exhibits a rich phase diagram, including a second-order transition to a 
phase where superconductivity and magnetism coexist. Results are compared with experiments on 
Cerium-based heavy-fermion superconductors. To study cold atom mixtures, we also simulate the 
BEG and two-component BEG models with a trapping potential. In the BEG model with a trap, 
there is no longer a first order transition to a true phase-separated regime, but a crossover to a kind 
of phase-separated region. The relation with imbalanced fermi-mixtures is discussed. We present 
the phase diagram of the two-component BEG model with a trap, which can describe boson-boson 
mixtures of cold atoms. Although there are no experimental results yet for the latter, we hope that 
our predictions could help to stimulate future experiments in this direction. 



I. INTRODUCTION 

Mixtures of 3 He and 4 He atoms exhibit a rich phase 
diagram, where besides a normal phase, there is a phase 
where 4 He is superfluid, and a phase separated region of 
superfluid 4 He and normal 3 Hei 1 - In 1971, Blume, Emery 
and Griffiths 2 proposed a model to describe such mix- 
tures. They simplified the continuous phase of the su- 
perfluid order parameter such that it could acquire only 
two values. Although they made this very rough ap- 
proximation and modelled the uniform system in a lat- 
tice, their results are very interesting. Qualitatively, 
they reproduced the right phases and the right orders 
of the phase transitions. Furthermore, if disorder is in- 
troduced by placing the mixture into aerogel, after some 
modifications^ the model can still yield the experimen- 
tally observed phase diagram. 4 

Here, we generalize this model to a two-component 
case in order to describe systems with two order param- 
eters and study the problem numerically, using Monte 
Carlo simulations. The motivation for the model we are 
proposing is twofold. Firstly, we would like to study con- 
densed matter materials like heavy fermions, high-T c su- 
perconductors, and organic superconductors. In particu- 
lar, we want to study the interplay between magnetic and 
superconducting ordering in these materials. Both order 
parameters are modelled as an Ising spin variable. Con- 
cerning the magnetism, we consider the ferro- and the 
antiferro-magnetic cases, and investigate also the effect of 
an additional magnetic field. We find that in the absence 
of a magnetic field, in the region where the two orders 
coexist, the system is always phase separated. When 
we add a magnetic field, we also find regions with mi- 
croscopic coexistence of the two phases. Secondly, we 
want to study mixtures of cold atoms. Cold atoms have 
emerged in recent years as an ideal simulator of con- 
densed matter systems. Because experiments with cold 
atoms are often carried out in a trap, we add a trapping 



potential to the model. This fact qualitatively changes 
the physics of the problem. For the case of a single com- 
ponent BEG model in a trap, the results are compared 
with experimental and theoretical work on imbalanced 
Fermi mixtures. For the case of the two-component BEG 
model, we make predictions for the phase diagram of 
boson-boson mixtures. 

The outline of this paper is the following: in section HT1 
we introduce the two-component BEG model, and inves- 
tigate it in the presence and absence of an external mag- 
netic field. The effect of a trapping potential is described 
in section [TTT1 In section HVl we compare the results with 
magnetic superconductors and cold atom systems. Our 
conclusions are presented in section fVl 



II. THE TWO-COMPONENT 
BLUME-EMERY-GRIFFITHS MODEL 

The BEG model was originally proposed to de- 
scribe superfluidity? 2 . The phase diagram found by 
Monte Carlo simulations exhibits large similarities with 
the phase diagram of 3 He- 4 He mixtures measured by 
experimentalists^ The main idea of studying superflu- 
idity with the BEG model relies on the U(l) symmetry- 
breaking of the ground-state wave function. For super- 
conductivity and Bose-Einstein condensation we have the 
same symmetry breaking, hence we can try to model 
these phenomena in the same way. 

Several physical systems exhibit two unequal symme- 
try broken phases simultaneously. A general Hamiltonian 
describing this class of systems reads 

Ti = -Ji ^2 CiVj - Ji X] s%s i + a i + cr " 

<ij> <ij> i i 

(1) 

where (cr^Sj) can take the values (0, 1), (0, — 1), (1, 0), 
and (—1, 0). This choice implies that only one kind of bo- 



2 



son can occupy each lattice site. D is an anisotropy field 
that controls the number of lattice sites with nonzero Oi. 
H plays the role of an external magnetic field, which may 
couple only to the order parameter describing a magnetic 
transition. The Hamiltonian ([1]) is appropriate for de- 
scribing phase transitions which require two order param- 
eters, one describing the ordering of the fraction of the 
system with nonzero a, the other one of nonzero s. This 
yields several possibilities, both fractions can model su- 
perfluidity, superconductivity, or (anti)fcrromagnetism. 
Possible applications could be magnetic superconductors, 
or two-component Bose-Einstein condensates. 

From now on, we will consider the fraction with non- 
negative a as describing magnetism, and s superconduc- 
tivity (preformed bosons that can Bose-Einstein conden- 
sate). Thus, <7, represents the spin of particle i and Sj 
the discretized phase of the wavefunction. Therefore, J\ 
can be both positive (ferromagnetism) and negative (an- 
tiferromagnetism) , but J2 has to be positive. We define 
the concentration, the ferromagnetic, antiferromagnetic, 
and superconducting order parameters as 

= -y yy» (3) 

i 

"»af l0 - = Jj ^(-l)Vi, (4) 

i 

m s = Jf^2 s i- ( 5 ) 

i 

Note that mf mjCT and m a f j(T can reach a maximum value 
of c, and m s of 1 — c. We define the ratio between the 
two coupling constants J2 and J\ as 

A. The Method 

We investigate this model by Monte Carlo simulations. 
To determine the location of second-order phase transi- 
tions, we performed simulations at constant concentra- 
tion, in which the elementary moves were flips of Sj and <7j 
or nonlocal spin exchanges. The location of the transition 
is then obtained from the peak location of the magnetic 
susceptibility. The locations of first-order phase transi- 
tions are obtained from simulations at constant temper- 
ature, with as elementary moves local flips of Sj and Cj, 
as well as same-site replacements of by <Xj and vice 
versa. A jump in the concentration c as a function of 
the anisotropy field D is then the signature of the phase 
transition. 

All simulations are performed on lattices with approx- 
imately 40 x 40 sites. Per point in the phase diagram, 
simulations were run over 3 ■ 10 5 to 3 • 10 7 Monte Carlo 
steps per site, depending on the correlation times. 



B. Zero magnetic field, H = 

In the absence of a magnetic field, the Hamiltonian ((T|) 
has ferro-antiferromagnetic symmetry. 

First, we consider K = 1. In this case, J\ = J2, and 
the shape of the phase diagram must be symmetric under 
the transformation c — » 1 — c. The results of the simula- 
tions are plotted in Fig. [T] We see that it indeed obeys 
this symmetry and exhibits four phases: a superconduct- 
ing phase (S), where the order parameter m s is nonzero, 
a ferromagnetic phase (FM), where mf miCr is nonzero, a 
phase-separated regime (PS) where the spins and the an- 
gular phases have formed ordered clusters, and finally 
the normal phase (N), in which there is neither order 
nor phase separation. Analogous to the BEG model, 
the transition from the phase-separated regions to other 
phases are first-order (dashed line), the other ones are 
second-order (continuous line). 
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FIG. 1: (color online) Phase diagram, temperature (in units 
of J\/ks) versus concentration, in the absence of a magnetic 
field. N indicates the normal phase, S superconductivity, FM 
ferromagnetism, and PS phase separation. Solid lines repre- 
sent second order phase transitions, dashed lines first order 
ones. Lines are guides to the eye. Snapshots of the simulation 
are shown. Black (white) represents o% = \ ( — 1), red (blue) 
Si = 1 (-1). 

Second, we consider the case K = 0.1. The results 
of the simulations are plotted in Fig. [3] We can under- 
stand the results as follows: J2 is much smaller than Ji, 
hence the spins will not pay much attention to the angu- 
lar phases, and the part of the phase diagram concerning 
the spins will be very similar to the BEG model. Be- 
cause J2 is so small, the phases will only order at very 
low temperatures (at zero concentration, the tempera- 
ture is ten times lower than the one at which the spins 
order at a concentration of one). If the concentration 
is slightly raised from zero, the system is already in the 
phase separated regime. All the states with a nonzero 
phase have clustered, and are not diluted by states with 
nonzero spin. Therefore, the critical temperature in the 
phase separated region will approximately remain con- 
stant. Because the temperature at which the angular 
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FIG. 2: (color online) Phase diagram in the absence of a 
magnetic field, for a relative coupling constant of K = 0.1. 

phases order is lower than the temperature at which 
phase separation begins, there is a phase separated re- 
gion in which the angular phases of the wavefunction are 
not ordered, which may appear unexpected at first sight. 
The transition within the phase separated regime, from 
the region where the angular phases are not ordered to 
the phase where they are ordered (superconductivity), 
is second-order. This is expected, because in the phase 
separated regime, all the phases have clustered, and the 
transition will be comparable with the transition in the 
Ising model, which is also second-order. 



C. Adding a magnetic field: the antiferromagnetic 

case 

If we apply a nonnegative uniform magnetic field to 
the system, the ferro-aiitiferromagnetic symmetry is bro- 
ken. We choose to consider the antiferromagnetic case 
here, because then there are two competing effects, the 
magnetic field tends to align the spins, whereas the ex- 
change interaction wants to order the spins antiferromag- 
netically. The magnetic field H will be measured in units 
of Ji. 

Kimel et al£- have studied the antiferromagnetic BEG 
model in the presence of a magnetic field, using Monte 
Carlo simulations. Their results at zero temperature 
suggest that the behavior of the system should be sep- 
arated into three qualitatively distinct regions, namely 
He [0,2], J? S [2,4] and H e [4,oo]. We consider here 
the cases K = 1 and K — 0.1 for values of H within each 
of these intervals. 



1. H=1.5 

First, we considered a magnetic field in the interval 
[0,2], namely H = 1.5. Both for K = 1 and K = 0.1, 
the results (not shown) are qualitatively the same as in 



the case of H = 0. This behavior was expected from 
the phase diagram of the single-component BEG model 
at zero temperature. Because the magnetic field tries to 
align the spins, the antiferromagnetic transition temper- 
ature is lower than in the absence of a magnetic field. 



2. H=2.5 

In the usual BEG model, the first-order phase tran- 
sition disappears in the presence of a magnetic field 
H G [2,4]. At zero temperature, there is a second-order 
phase transition between a state with oi = at every 
site, and a checkerboard phase, where one sublattice has 
<ji = at every site, and the other one <7j = — 1. There 
is also a transition between the checkerboard state, and 
an antiferromagnetic phase, but this transition is absent 
at nonzero temperature^ 

For K = 1, the behavior of the two-component BEG 
model is still very similar to the case H = 0. For K = 0.1, 
the first-order phase transition disappears, and therefore 
there is no phase-separated region, see Fig. [3] We do ob- 
serve an antiferromagnetic and a superconducting phase, 
but it is not clear from the figure whether the two phases 
overlap. To better understand this low-T intermediate 
regime, we also simulated the problem at a relative cou- 
pling strength of K — 0.5. In Fig. [4] we clearly observe 
that there is a region where antiferromagnetism and su- 
perconductivity coexist, without true phase separation, 
since the first-order phase transition has disappeared. 
What is also interesting is that at zero temperature this 
region begins at a nonzero concentration, and ends at 
a concentration smaller than one. When there is phase 
separation, this coexistence region always begins at c = 
and ends at c = 1. 

H=2.5 K=0.1 
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FIG. 3: Phase diagram at a magnetic field H — 2.5 and 
relative exchange strength K — 0.1. N denotes the normal 
phase, S superconductivity and AF antiferromagnetism. 
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FIG. 4: (color online) Phase diagram at a magnetic field H = 
2.5 and relative exchange strength K — 0.5. There is a region 
where superconductivity and antiferromagnetism coexist, but 
where there is no true phase separation. 

3. H=5 

In the original BEG model, when the magnetic field is 
increased to a value higher than H = 4 at zero tempera- 
ture, antiferromagnetism totally disappears because the 
spins tend to align with the magnetic field£ The system 
is therefore magnetized, but not because of the nearest- 
neighbor interactions. Therefore, this is not really ferro- 
magnetism, but for the sake of simplicity, we denote it 
like this. For the case of K = 1, we observe a phase with 
ferromagnetic and superconducting ordering, and a fer- 
romagnetic phase (not shown). For K = 0.1, we find an- 
other interesting phase, namely a ferromagnetic checker- 
board phase, consisting of two sublattices, see Fig. \5\ At 
the first sublattice, all sites are randomly occupied by 
phases with a value of s, = 1 or Sj = — 1. At the second 
one, all sites are occupied by the spin that is favored by 
the magnetic field, cr^ = — 1. This phase is most likely to 
occur at a concentration of c = 0.5 because in this case 
a perfect checkerboard is possible. 
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FIG. 5: (color online) Phase diagram at H — 5 and K = 0.1. 
FM denotes ferromagnetism, S superconductivity and CFM 
denotes checkerboard ferromagnetism. 



III. ADDING A TRAP POTENTIAL 

A. The Blume-Emery-Griffiths model 

Because experiments with cold atoms are often car- 
ried out in a trap, we will add a harmonic potential to 
the original BEG Hamiltonian, to describe mixtures of 
fcrmions and bosons in a trap. In general, the potential 
felt by the bosons is different from the one felt by the 
fermions, what implies that we must include two terms, 

+ M + «/ X>? + % 2 )(i - (7) 

i i 

Here, Xj and y - L are the horizontal and vertical distances 
of site i, measured from the center of the lattice, in lattice 
units, and a& and a / measure how much the bosons (the 
states with = ±1), and the fermions (the states with 
a, = 0) feel the influence of the trap. If a& = a/, this 
term is constant, and the phase diagram is not modified. 
We will consider the case a& > aj, which is the most 
relevant experimentally. Using the hard core constraint 
of + s 2 = 1 , we can then rewrite this term and add it to 
the BEG Hamiltonian, thus obtaining 

<ij> i i 

where a = aj, — df. This means that, effectively, the 
bosons will feel a stronger tendency to go to the center 
of the trap. 

In the limit of a — > oo, all the states with <7j = ±1 will 
cluster in the center of the trap, and therefore the order- 
ing temperature will be the same as in the Ising model. 
Note that the maximum value of the extra term in the 
Hamiltonian will depend on the size of the lattice. This 
way of including the trapping potential is comparable 
with the work of Gygi et al.^ where a spatial-dependent 
chemical potential was added to the Bose-Hubbard model 
in order to describe bosonic atoms in an optical lattice. 

We simulated the new model using the same proce- 
dures as for the BEG model and the two-component 
BEG model. The results for three different strengths 
of the trapping potential are plotted in Fig. [51 In the 
BEG model without a trap, there is a second-order phase 
transition from a normal state to an ordered state, and a 
first-order phase transition to a phase-separated region 2 . 
For the three values of a considered here, we do not find 
a first-order phase transition any more. A part of the 
first-order phase transition line disappears, and a part 
changes into a second-order one. 

We see that for a small difference between the trap po- 
tential felt by the bosons and the fermions, a/J = 0.001, 
the transition temperatures are very similar to the case 
without a trap. For a large difference, a/J = 0.1, the 
transition temperatures approach the transition temper- 
ature of the Ising model for almost all concentrations, as 
expected. When the states with Uj = ±1 are ordered, we 
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FIG. 6: (color online) Phase diagrams of the BEG model with 
a trapping potential. N denotes the normal, unordered state, 
C the condensed phase, in which the sites with <?i = ±1 arc 
ordered. Lines are guides to the eyes. Snapshots are shown, 
where blue represents <Xi = 1, red cti = —1, and white at — 0. 



lattice of size L = 41 will start to feel the potential if 
k B T/J ~ 800a. For a/J = 0.1 and a/J = 0.01 this 
results in ksT/J ~ 80 and ksT/J ~ 8, respectively, in 
both cases much higher than the temperatures we are 
interested in, because ordering starts around ksT / J « 
2.4. Therefore, the single particles will experience the 
influence of the trap in the entire temperature range of 
Fig. [5] (b) and (c). For a/J — 0.001, a single particle 
will feel the trapping potential for temperatures lower 
than kBT/J ~ 0.8. However, for higher temperatures the 
system already orders, and therefore there are some large 
clusters that according to Eq. @ will feel the potential 
already at much higher temperatures. This reasoning is 
in agreement with the snapshots in Fig. 0(a). For a/J — 
0.001, we clearly observe the influence of the trap when 
the states with <t^ = ±1 have clustered. In the disordered 
state, the influence is less visible. For a/J = 0.1 and 
a/J = 0.01, we indeed see the influence of the trap for 
all temperatures, even in the disordered state. 



B. The two-component Blume-Emery-Griffiths 
model 



Analogous to the previous subsection, we will also add 
a trapping potential to the two-component BEG model. 
In the latter, both the states with <7j = ±1 and Si = ±1 
describe bosons, that both can condense. Therefore, this 
model can be applied to study cold atoms mixtures with 
two species of bosons. We will consider the realistic case 
that the two species feel different trapping potentials. 
Therefore, we add the extra terms 



2\ 2 



a s ^2(Xi + yf)s 2 



(10) 



to the Hamiltonian. Because at every lattice site <j| + 
s 2 = 1, we can rewrite this term and add it to the two- 
component BEG Hamiltonian, to get 

Ti = -Ji ^2 <T t (J j ~ h 53 SiS i + D y^. a i 

<ij> <ij> i 

+aY J {.x 2 l +y 2 l )al (11) 



will not speak of a superfluid state, but of a condensed 
state, because we now consider bosons in general. 

It is important to estimate at which temperature the 
system starts to feel the influence of the trapping po- 
tential. Let us assume that a cluster of size m feels the 
potential when the energy difference between the state 
with this cluster in the center and in the corner of the 
lattice is of the order ksT/J. For a lattice of size L 2 , 
this estimation results in 



maL 
2J 



k R T 



(9) 



In this approximation, a single particle (m = 1) in a 



where a = a„ — a s is now the difference between the 
potentials felt by the two species of bosons. Now, the 
bosons with Cj = ±1 have a stronger tendency to go to 
the center of the lattice. 

The results of our simulations are plotted in Figs. [7] 
and [8] We considered two different strengths of the trap- 
ping potential, and two different ratios of the coupling 
strengths of the bosons, namely K — J2/J1 = 1 and 
K = 0.1. For K = 0.1, the right part of the first- 
order phase transition disappears and the left one be- 
comes second-order (see Fig. [7]), whereas for K = 1 both 
left and right parts of the first-order phase transition are 
converted into second-order (see Fig.[H]). 
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FIG. 7: (color online) Phase diagrams of the two-component 
BEG model with a trapping potential. TV denotes the normal, 
unordered state, C a and C a the phases where the bosons rep- 
resented by the state with <7j = ±1, respectively Sj = ±1 
are condensed. Lines are guides to the eyes. Snapshots are 
shown, where black and white represent <Xi = ±1, and red and 
blue Si — ±1. 



K=1 



with Si = ±1 at the corners. Therefore, for all concentra- 
tions, the system behaves as two uncoupled Ising models. 
In the case of K = 1, we see indeed that the transition 
temperatures for both species approach the Ising transi- 
tion temperature. For K = 0.1, because Ji is ten times 
smaller than J\ , one of the species will order at the Ising 
transition temperature, and the other one at one tenth 
of the Ising transition temperature. 

To find the temperature at which the system starts 
to feel the presence of the trap, we can make the same 
analysis as in subsection IIII Al Also here, we see in the 
snapshots of Figs.[7]and[5]that for a/ J\ = 0.1 (not shown) 
and a/Ji = 0.01, the system always feels the influence of 
the trap, and for a/Ji = 0.001, it does only when the 
system is ordered. If we inspect Fig. [8] (a), we see that 
there is a phase C s in which the bosons represented by 
,s; = ±1 are ordered, but the bosons represented by Oi = 
±1 are not. This is somewhat surprising. A reason for 
the occurrence of this phase is that when all the bosons 
that have the tendency to go to the center of the trap 
have clustered there, automatically also the other bosons 
have clustered at the edge. Therefore, they can have 
nearest-neighbor interactions, and they can easily order. 
It remains to see whether such a phase indeed occurs in 
experiments. From the theoretical point of view, it would 
be interesting to also allow for states with <Ji = Si = 0, 
to verify the stability of this phase, when we relax the 
constraint that every lattice site must be occupied by one 
of the bosons. Note that for small enough concentrations, 
this phase will always occur, since the bosons s are hardly 
diluted by the bosons <j. 
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IV. COMPARISON WITH EXPERIMENTS 



A. Magnetic superconductors 
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FIG. 8: (color online) Phase diagrams of the two-component 
BEG model with a trapping potential. The notation used is 
the same as in Fig. [7] 



FIG. 9: Phase diagram of CeCo(Ini_ I Cd I )5. The figure is 
extracted from Ref. M. 



There are several examples of Cerium-based supercon- 
In the limit of a — > oo, all the sites with Oi = ±1 will ductors, for example CeCoIns and Celrln5, as well as an- 
have clustered in the center of the lattice, and all sites tiferromagnets that contain this element, like CeRhlns, 
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CeCoCds, CeRhCd5, and CeIrCd 5 . Let us consider 
CeCoIn 5 and CeCoCds. These two materials have two 
elements in common, Ce and Co, and differ in the third 
element. By doping CeCoIns with Cd on the In site, 
we can change the superconductor CeCoIn.5 into an an- 
tiferromagnet. There are more of these Cerium-based 
pairs, and therefore, this class of materials is appropri- 
ate for studying the interplay between superconductivity 
and magnetism. 

Let us consider the heavy fermion superconductor 
CeCoIns, with Cadmium doping on the In-site. This ma- 
terial has the highest superconducting transition temper- 
ature (T c — 2.3K) of all heavy fermions, and its electronic 
structure is quasi-2Di£ Nicklas et al& and Pham et aL— 
determined the antiferromagnetic and superconducting 
onset temperatures of this material as a function of dop- 
ing by elastic neutron scattering, specific heat, and resis- 
tivity measurements. Their results are plotted in Fig. [5] 
For experimental details we refer the reader to Ref. 8. 
The phase diagram of CeCo^ni-^Cd^s shows three or- 
dered phases: a superconducting phase, a commensurate 
antiferromagnetic phase, and a region where supercon- 
ductivity and antiferromagnetism microscopically coex- 
ist. 



is pressure. Park et alr^ determined the phase diagram 
of CeRhlns in temperature-pressure space with and with- 
out a magnetic field. Without a magnetic field, they also 
found this abrupt disappearance of the incommensurate 
antiferromagnetic order at S±. However, when they ap- 
plied a field of 33 KOe, the line of the magnetic ordering 
temperature went smoothly down to zero at 62- Such a 
phase diagram shows many similarities with Fig. [5] if we 
identify pressure with inverse concentration in our model. 
Indeed, for an external magnetic field of H — 2.5 and a 
relative exchange strength K — 0.5 (see Fig.|4|), the phase 
diagram shows the same three ordered phases. Further, 
the coexisting phase is not phase separated. 

Finally, we consider the compound CeIr(Ini_. c Cd a ;)5, 
see Fig. Eland Ref.©. For this material, it is not clear if 
there is a region where superconductivity and magnetism 
coexist. If there is such a region, it is in a small doping 
interval. The phase diagram of this material strongly re- 
sembles the phase diagram of the two-component BEG 
model with an external magnetic field of H = 2.5 and 
a relative coupling strength of K — 0.1, see Fig. [3] Al- 
though this experiment was also carried out without an 
external magnetic field, we only find similarities with our 
model in the presence of a magnetic field. 
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FIG. 10: Schematic phase diagram of unconventional super- 
conductors in temperature-control parameter space. AF de- 
notes antiferromagnetism, S superconductivity and NFL a 
non-Fermi liquid. Experimentally, antiferromagnetism often 
disappears abruptly at some critical value <5i of the control 
parameter, although one would expect a magnetic quantum 
critical point at some value S2 of the control parameter. 



It is interesting to observe that in this material antifer- 
romagnetism suddenly disappears at the point where the 
onset temperatures for superconductivity and antiferro- 
magnetism are equal. This feature, however, may change 
in the presence of an applied magnetic field. In Fig. [TU] 
we see a schematic phase diagram of unconventional su- 
perconductors, in temperature-control parameter space. 
In the case of CeCo(Ini_ a; Cd a ;)5, the control parameter 
would be doping. Another example of such a parameter 
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FIG. 11: Phase diagram of the heavy fermion 
Celr(lni_a;Cd:r)5. The figure is extracted from Ref. 0. 



B. Cold atom systems 

In 2006, two experimental groups, the group of 
Ketterle at MLTjii and the group of Hulet at Rice 
University^ have performed experiments with imbal- 
anced ultracold 6 Li atoms in a trap, and obtained con- 
tradictory results. The MIT group measured a transition 
between a normal and a superfluid phase at a polariza- 
tion of P k 0.70, whereas the group at Rice University 
observed a transition between two superfluid phases at 
P w 0.09. Here, P measures the imbalance between the 
spin-up and the spin-down atoms, 



iV T + Ni 



(12) 



Gubbels et ali^ have set up a theoretical model to de- 
scribe these imbalanced Fermi mixtures and determined a 
general phase diagram in temperature-polarization space 
that can explain the observations of both groups. The 
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topology of their phase diagram shows large similarities 
with the phase diagram of the BEG model. We can un- 
derstand this resemblance as follows. In the BEG model, 
the concentration c is the fraction of lattice sites with 
<Ji = 0, and thus the fraction of the system that cannot 
condense. The polarization P is a measure for the differ- 
ence of the atoms in the spin-up and the spin-down state, 
and thus for the number of fermions that remain after 
the others have paired. The atoms with spin up and spin 
down will form pairs, and such a pair can be described 
as a boson. Therefore, the polarization is also a measure 
for the fraction of the system that cannot condense, and 
the concentration can be mapped onto the polarization. 
We can identify the paired atoms, the preformed bosons, 
with the states Oi — ±1, and the remaining fermions with 
a* = 0, see Fig. [£2] 



2.5 




FIG. 12: Phase diagram of the original BEG model, obtained 
by Monte Carlo simulations. N denotes the normal phase, 
SF superfluidity. Lines are guides to the eyes. The transition 
between the normal and the superfluid state is second-order, 
the transition to the phase separated regime is first-order. 

The experiments with 6 Li are carried out in a trap, 
and the theoretical model of Gubbels et al. only includes 
the presence of the trap by using the local density ap- 
proximation. Now, we would like to compare their phase 
diagram with our results of the BEG model with a trap- 
ping potential, in Figs. [7] and [5] Although in the case of 
imbalanced fermions the frequency u> of the optical trap 
felt by the pairs of fermions (bosons) and the remaining 
unpaired fermions is the same, the mass of the bosons 
is twice as large, and the potential constant at is thus 
larger than af. This means that the comparison must 
be made with the BEG model in a trap. In this model, 
the first-order phase transition, measured by a jump in 
the concentration as a function of the anisotropy field 
D has disappeared, thus there is no true transition to 
a phase-separated regime. However, if we inspect the 
snapshots, we see that for low enough temperatures, or 
large enough trapping potential, there still is a clear sep- 
aration between the condensed bosons and the fermions, 
suggesting some kind of effective phase separation. We 



note that in experiments, phase separation is measured 
by inspecting the radii of the clouds of the atoms in the 
different hyperfine states, and not by a jump in some or- 
der parameter^ Our results thus suggest that the mea- 
sured different radii are not per se an evidence of a true 
thermodynamic phase separation. Further experiments 
are required to clarify this issue. 

Although our model describes qualitatively the exper- 
imentally observed phases, it cannot capture the fine de- 
tails of recent experimental results. Studies by Shin et 
al^ indicate that there is no superfluid phase, or phase- 
separated phase for polarizations above P 0.36. By a 
quantum Monte Carlo approach, Lobo et al^ predict a 
phase transition between a normal and a superfluid state 
at a polarization of P ~ 0.39 at zero temperature, and 
Gubbels and Stoof 1 — recovered this results using a Wilso- 
nian renormalization group theory. 



V. CONCLUSIONS 

We simulated a two-component extension of the BEG 
model without an external magnetic field and determined 
the phase diagram in the concentration-temperature 
space. In the region where magnetism and superconduc- 
tivity coexist, the system is always phase separated. We 
added a magnetic field to our model, and considered the 
antifcrromagnetic case. In this case, we also find phase 
diagrams with true coexistence of two ordered phases. 
These diagrams are comparable with the phase diagram 
of doped heavy fermions in the presence of a magnetic 
field. 

In order to describe cold atom systems, we added a 
trapping potential to the BEG model, and our extension 
of this model. The added potential changes the phase 
separation regime conceptually. We cannot speak any- 
more about true phase separation, but more about a 
crossover to a phase separated region. We argue that 
the BEG model with a trapping potential can be used to 
model imbalanced Fermi mixtures. However, there are 
still quantitative differences with experiments, which our 
model is not able to cover. We also made predictions 
for the phase diagram of boson-boson mixtures based on 
our simulations of the two-component BEG model with 
a trapping potential. Although there is no available ex- 
perimental data on boson-boson mixtures, we hope that 
our work can motivate further studies in this direction. 
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